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ABSTRACT 

We discuss various mechanisms that can lead to interband sign reversal of the order parameter in a multi- 
band superconductor. In particular, we generalize Abrikosov-Gor'kov solution of the problem of weakly coupled 
superconductor with magnetic and nonmagnetic impurities on the case of arbitary order parameter anisotropy, 
including extreme cases as d— pairing or interband sign reversal of the order parameter, and show that interband 
scattering by magnetic impurities can stabilize an interband sign-reversal state. We discuss a possibility of such 
state in YBa2Cu307 in the context of various experiments: Josephson tunneling, neutron scattering, isotope effect 
measurements. 
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1. In the last years new experiments with high-temperature superconductors revived the interest to the strongly 
anisotropic superconductivity. Under this name we mean superconductors where the order parameter A is strongly 
different on different parts of the Fermi surface. It includes the cases when strong anisotropy of A is implied by 
symmetry of the superconducting state, like d-wave pairing, one-band superconductors with a strong anisotropy, 
but the same symmetry as in the normal state (so-called extended s pairing), or multiband superconductors with 
substantially different order parameters in different bands. It does not include superconductors with anisotropic 
electronic structurewhich does not result in anisotropic order parameter. 



In this paper we shall concentrate on the case of multiband superconductors, more precisely on the case of 
extreme interband anisotropy when the sign of the order parameter is different in different bands. We shall discuss 
various conditions which can lead to such a state. First, we shall address the question whether a combination of 
phonon-induced attraction and homogeneous Coulomb repulsion may lead to interband sign reversal of the order 
parameter. Second, we shall derive a generalization of the Abrikosov-Gor'kov formula for superconductivity in 
dirty superconductors!!] for arbitrary anisotropy of the order parameter, and will show that impurity scattering 
does not necessarily reduce the anisotropy of the superconducting state. Moreover, we shall point out cases when 



impurity scattering can stabilize a solution with an interband sign reversal. Finally, we will describe a model when 
pairing appears due to electron-paramagnon scattering and the resulting state has an interband sign reversal of 
the order parameter. Finally, we will discuss various experiments in YBa2Cu307, and show how they can be 
explained with the interband sign reversal concept. 

2. The extension of the BCS theory for two or more superconducting bands was first worked out by SuhL 
Matthias, and Walkercl and independently by Moskalenko,0 and later elaborated on by many. It was realizeoa 
that the fact that several bands cross the Fermi level is not sufficient to have considerable many-band effects in 
superconductivity. Only when the bands in question have a very different physical origin, can a substantial effect 
appear. It has also been realized that the anisotropy, if any, appears already in the weak coupling (BCS) solution, 
while the proper strong coupling treatment addsno qualitative features. 

Many high-T c cuprates not only have multisheet Fermi surfaces, but the actual bandsat the Fermi level are 
qualitatively different. In particular, YBa2Cu307 is known to have four sheets of the Fermi surface, all four 
havingdiffcrcnt physical originaJ: One is formed by the chain pda states (seen by positron annihilation), another 
is an apical oxygen band (seen in de Haas- van Alphen experiments), and the last two are bonding and antibonding 
combinations of the two pda plane bands (seen by angular-resolved photoemission) . Basing on the richness of 
the band structure of YBa2Cu30y, several groups pointed out that at least the two-bandj3 or probably the whole 
four-bandJ3'EI picture should be used to describe superconductivity in this system. Various experiments have been 
interpreted as indicating two or more different superconducting gaps. 



We shall now remind the basic equations of the multiband BCS theoryaH: The Hamiltonian has the following 
form: 
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where e,^ is the kinetic energy in the i-th band, c* ka and Ci^a are corresponding creation and annihilation 
operator, and gtj is the averaged pairing potential. 

The order parameter A on the i-th sheet of the Fermi surface is given by the equation 
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if the cut-off frequency ujd is assumed to be the same for all sheets. T c is defined in the usual way by the effective 
coupling constant, log(27*w d/kT c ) — 1/A e //, 7* ~ 1.78. The effective coupling constant A e // in this case is 
simply the maximal eigenvalue A max of the matrix Ay = gtjNj, where Nj is the density of states at the Fermi 
level (per spin) in the j-th band. Ay plays the role of the coupling constant A in the one-band BCS theory. Note 
that conventional (isotropic) A is also defined in terms of Ay: A = J2ij AyA^/A = XiNi/N, where the mass 
renormalization for the i-th band, as measured, for instance, in de Haas- van Alphen experiments, is Xi = Y]. Ay, 
and N = yV Nj. Obviously A e // > A, which means that due to larger variational freedom T c in the multiband 
theory is always larger than in the one-band theory. The two are equal in isotropic case, i.e. when gy does not 
depend on An instructive example of the opposite case is the two-band model with An = A22 = A > 0, 
A12 = A21 = —A. Then A = 0, while A e // = 2A. Note that the last value is the same as when A12 = A21 = A. 
The physical reason is that although there is no solution of Eq.[l] with Ai > 0, A2 > 0, there is an obvious solution 
with Ai = — A2 7^ 0. Near T c , the solution of Eqjl] is A2/A1 = (A e // — An)/Ai2, demonstrating directly that 
the sign reversal of the order parameter, A2/A1, takes place when nondiagonal matrix elements A12 and A21 are 
negative. One can easily check that Eq.|l| may have a superconducting solution even for all c/y < 0, i.e., when no 
attractive interaction is present in the system. The condition for that is \gi2\ > (\gi\N 2 + 1 A| )/2A r 1 7V 2 . This is 
similar to the well-known fact that in a system with repulsion the superconductivity with higher angular momenta 
(p, d) is possible, because of the sign reversal of the order parameter. The main difference is that in the example 
above the symmetry of the superconducting state is the same as of the normal state. Below we shall demonstrate 
that even a fully attractive interaction gij > can lead to the sign reversal if (a) interband pairing interaction is 
weaker than Coulomb pseudopotential, or (b) there is strong interband scattering by magnetic impurities. 
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Figure 1: Critical temperature and superconducting gaps in a system with induced superconductivity in the 
second band, as a function of Coulmb pseudopotential. 



If g's are electron-phonon pairing potentials, then Eq.|l| should be corrected for a Coulomb repulsion, which 



can be readily doneo by substituting g^j 



9ij 



U* 



U* , where the effective Coulomb repulsion U* 



is logarithmically renormalized in the same way as in one-band superconductivity theory (U* is assumed to be 
independent on A direct consequence of that is that if the interband electron-phonon coupling is weak, 

the situation with a negative gap, gij — U* < 0, can easily be realized because of the interband repulsion. We 
illustrate that by numerical calculations presented in Fig{j]. In these calculations the following parameters had 
been used: <?i2 = 522 = 0, Ni = 4JV2, and 511 = N± so that to have A = 1. 



Several facts draw attention: First, in this model T c decreases with the increase of /1* = U*N substantially 
slower than in a one-band case when T c — > when fj,* — > A. Second, the order parameter induced in the second 
band (not superconducting by itself) is always negative; Its absolute value reaches maximum when |Ai| = IA2I, 
i.e., at U* = g pp Np/2(Np — N c ). Further below we shall discuss the relevance of this situation to YBa2Cu30y. 



3. Now we shall analyze the influence of impurity scattering on the sign and anisotropy of the order parameter. 



Following the standard way of including the impurity scattering in the BCS theory ,§ one writes the equations for 
the renormalized frequency uj n and order parameter A n (n is the Matsubara index), which completely define the 
superconductive properties of the system: 

A Jn = Aj + Y^^hiJ'-lh') (2) 
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The only difference from a specific case of multiband superconductivity (i.e., interband anisotropy), considered 
in a number of papers (see, e.g.— .Refs) is that instead of the band indices in || appear indices J, J', labeling 
Allen's Fermi surface harmonics. E£l Note that for a spherical (cylindrical) Fermi surface Allen's harmonics reduce 
to spherical harmonics and J to {L,m} or just m. For the interbancLanisotropy without intraband anisotropy J 
can be chosen to be the band index (so-called disjoint representation!!!]). Other notations in EqB have their usual 



meaning: uj n = (2n + 1)ttT, Qj n — y £>j„ + Aj n , Ijj 1 is the scattering rate matrix bynonmagnetic impurities, 
and Xjj, is the same for magnetic impurities. Coupling matrix A is defined in the same way as Allen'smatrix 
Xjj,B^jj' = V^Njj,(0). 

Linearizing Eqs. ^ with respect to A (an approximation valid near T c ), we get: 

A' Jn (l + Gj/2w n ) =Aj + J2 *!j, n Tjj,/2ui n => A' Jn = £ Aj,{8jj, + gjJ ,/2uj n )-\ 



where we have defined 



Z Jn = l + Gj/2u n ; A Jn = Aj + J2 Aj ' n Tjj, 

A' /n = Aj n /Zj n ; Gj = ^7,/j'; Z Jn =uij n /w n 

j' 

Tjj' = 7j./' - Ijj'] gjj' = Sjj'Gj - Yjji. 



From this immediately follows that 
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which after a usual trick with subtracting the clean limit, g — 0, and extending summation to infinity (a useful 
matrix formula is (l+A) _1 =I-A(I+A) _1 ), gives 

Aj = J2 A JJ'[LSj'k ~ Xj>k]A k ; Xjj. = ^9kJ' /2 a)" 1 (Wn^/fc + gjk/ty' 1 , (3) 
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where L = ln(2"/* lj d / ftT c ) . By introducing the eigensystem of g, gjj' — ^2 k Rj^dkRkj' , we can express X in 
terms of the difference between the two incomplete gamma- functions: 

Xjj, =Y J R: J lY. bJ n ldk ^ + dkl2Y X R k j, =Y,R-jlx{d k )RkJ>i 
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with x( x ) — "0(1/2 + x/2ir) — -0(1/2), which is the standard definition of the matrix function X =x(g)- 



This result is analogous to the classical one of Abrikosov and Gor'kovJH but includes arbitrary anisotropy. 
Now solving (||) for L, we find: 

Aj = Y,^j 1 k+Xj k )- 1 LA k , (4) 
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which means that now T c is defined by the maximal eigenvalue of the effective matrix A. e tf = (A -1 + A) -1 . 
As can be seen immediately from Eqs.(||-^|) and the definition of g,jj', diagonal nonmagnetic scattering rates 

711 have dropped out from Eq.(^). This is the manifestation of Anderson theorem for an anisotropic case: such 
scattering does not influence T c (in a considered Born limit), as in isotropic superconductors. As will be discussed 
below, this argument works only for diagonal in J, J' non-magnetic scattering, while all other are, in principle, 
pair-breaking. 

In the second order in A (assuming that AX is small), 

A eff = A - AAA. (5) 

If we recall that A forms the eigenvector of A for the maximal eigenvalue A e //, we can immediately write the 
lowest-order correction to A e //: 

SX eff =-£ AjXjj, A j, /L 2 J2A 2 j (6) 
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It is illustrative to consider explicitly a two-band case for weak scattering (7^ , 7^ <C T c ) ; the effective matrix 
is then given by: 

A e// =A-— A.f 2 7"+^+Ti2 712-712 Va. (7 ) 

11 8T c0 V 7ii-72i 2 7 | 2 + 7^ + 721 ; w 

When all A's are equal (isotropic case), the standard Abrikosov- Gor'kov result is recovered: SX w — 7rA 2 (7f 1 + 

712 + 7li + 722)/87co. The main point of the AG theorytl is that 7 s enters equations for w and A with opposite 
signs. That is why the magnetic impurities appear to be pair-breakers, and the non-magnetic ones not. The 
above solution shows explicitly that in anisotropic case of Eqs.|| only intraband non-magnetic scattering does not 
influence T c (7^ drop out). 

An interesting special case is Ai2,A2i <C An,A22- Then in the effective A matrix the nondiagonal element 
h^J = A12 + vrAnA 2 2(7i2 — 7i2)/87 7 o becomes negative, if 7*3 —712 > 8T , c oAi 2 /AiiA 2 2. As discussed above, this 
situation will lead to sign reversal of the order parameter. One can say that, when attractive interband coupling 
is relatively weak and the magnetic interband scattering is strong, the system chooses to have two gaps of the 
opposite signs, losing in pairing energy, but avoiding the pair-breaking due to interband scattering. 

4. Maybe the most interesting example of an interband sign reversal of the order parameter is supercon- 
ductivity due to electron-paramagnon interaction, in other words, superconductivity due to dynamic exchange 
of spin fluctuations. For simplicity, we consider the case of singlet pairing. In this case the pairing potential 
is Vij, kl = J dRdli' J2 a f3-ys(i' a \J( r ~ ^■) <T a0\jP)x('^- — R')(^7l^( r _ R)vj5\lS) where J is exchange interaction 
and x = (S(R)S(R. )) is spin-spin correlation function. Contrary to the electron-phonon interaction, this inter- 
action is repulsive (for triplet pairing, it would be attractive). Thus it cannot lead to a superconducting solution 
where the order parameter has always the same sign. A possible solution is the well-known <i-wave state, where the 
order parameter changes sign upon rotation by 7r/2 in the momentumspace. Another possibility is to have order 
parameter of different sign in different bands. In this case the interaction must be small in the intraband channels 
and strong in the interband channel. It looks on the first glance to be an artificial, unnatural condition. However, 
in some cases it emerges quite naturally. A good example is a bilayer with spin fluctuations antiferromagnetically 



correlated between the layers. This system models some superconducting cuprates, like YBa2Cu30y_5. If the one- 
electron tunneling between the layers is larger than superconducting gap, two bands are formed and well defined: 
bonding (symmetric) and antibonding (antisymmetric). By symmetry, only electron states of different parity 
can interact via exchange of antiferromagnetic (thus antisymmetric) spin fluctuations. Spin-fluctuation induces 
superconductivity in YBa2Cu307_a with sign reversal of the order parameter was quantitatively investigated in 
Ref.jB 

5. Now we shall briefly discuss some experimental implications of the interband sign reversal of the order 
parameter, if it occurs in YBa2Cu3C>7. It is believed that four different bands cross the Fermi level in YBa2Cu30y: 
two plane bands, which are bonding (B) and antibonding (A) combinations of the individual planes' states, the 
chain (C) band, and a small pocket formed mainly by apical oxygen states (which is not discussed here). We 
adopt here the point of view that the superconductivity originates in the plane bands, and is induced in the chain 
band by interband proximity effect. Thus we shall consider a superconducting state which is characterized by 
the order parameters of the opposite signs in the bonding and antibonding bands. For the following discussion 
we shall also mention that, according to calculations,!] the chain band is very light, so that its contribution to 
the total density of states is small (~ 15%), while its contribution in the plasma frequency ojp y oc N(Q)v Fy is 
considerable (~ 50%). These finding are confirmed by the experiment: Maximal Fermi velocity was calculatcdLL3 
to be ~ 6 x 10 7 cm/s and corresponds to the point where the chain Fermi surface crosses the T— Y line. This value 
agrees wall with the Raman experiments. Ej Calculated plasma frequency anisotropy Wp^/w^ ~ 1.75, as discussed 
in Ref.,E£l is in agreement with the optical and transpprt measurements. Band A is, according to calculations, 
rather heavy, which also agrees with the experiment .E£E3 Band B is light again. Both A and B bands are nearly 
tetragonal. Their relative contribution to the normal-state transport is defined by the partial plasma frequencies 
(Table I). Importantly, bands A and C at q z = can cross by symmetry, for instance they are degenerate with 
e = Ep at q=(~ 0.87r/a, « 0.2tt/B, 0). For all q z ^ these bands hybridize. This is the reason for YBCO being 
the most three-dimensional of all high-T c cuprates. An extremal orbit in q z = plane, which appears because of 
the A — C hybridization, has been seen in de Haas-van Alphen experiments.EZI 

Which of the three possible mechanisms (if any) can be operative in YBa2Cu30y and produce interband 
sign reversal? It appears that any of the three can. First, it is possible that the phonon-induced interaction is 
stronger in the intraband than in the interband channel. For instance, one of the key feature of the electronic 
structure of YBa2Cu307 are extended saddle points at van Hove singularities. It-acas shown that extending of the 
saddle points is directly related to the geometry (warping) of the Cu02 planes x3 Saddle points appear only in 
the antibonding band. Thus, one may expect the saddle points regions to interact particularly strongly with the 
phonons of the even parity (gerade), with correspondingly weaker intraband interaction. Second, the main type of 
magnetic defects in YBa2Cu307 are localized magnetic moments on Cu, which appear near O vacancies and other 
imperfections. In view of very strong interplane antiferromagnetic correlations of the Cu spins in underdoped 
compounds, it seem plausible that each defect-related Cu moment in the fully doped system induced magnetic 
moment at the nearest Cu site in the next plane. Such pair of spins will have exactly that symmetry, that is 
needed to stabilize the superconductivity with the opposite signs of the order parameter. Third, the popular 
model of the spin-fluctuation superconductivity, which is usually associated with the d-pairin^i in fact always 
has another solution, namely that with the angular s-symmetry and interband sign reversalJlj Which of the 
two solutions is more stable depends on whether or not the spin fluctuations are antiferromagneticaly correlated 
between the planes. If they are, the second solution is always more stable. Correspondingly, if one accepts the 
results of the numerous neutron scattering experiments, which do show the interplane correlation, and assumes 
the spin-fluctuation induced superconductivity, the resulting state is one with the interband sign reversal. Note 
that this is true both for the overdamped magnetic excitation deduced by Pines and co-workers from the NMR 
data, and for the finite-energy paramagnon as deduced from the neutron scattering. 

Now we shall briefly discuss some experiments, probing the symmetry and/or anisotropy of the superconduct- 
ing state. It is instructive to single out those few experiments that deal with the symmetry of the order parameter. 
Those are Josephson tunneling experiments and neutean scattering experiments. Both were discussed in detail, 
in the context of the interband sign reversal, in Rcfs.1121 and,E2l respectively. Numerous Josephson tunneling ex- 
periments, performed in the last two years, show that the phases of the tunneling current along crystallographic 



a and b directions are shifted by it. A natural interpretation is in the form of the one-band d x 2_ y 2 state. As it 
was pointed out in Ref.Jlj because of the strong difference in effective masses, the chain band in YBa2Cu30y 
contribute substantially into the transport in b direction. Thus, the relative phase of the tunneling currents in 
the interband sign reversal model depends on whether the sign of the order parameter in the chain band is_the 
same as in the bonding band (which dominates in-plane transport) or as in the antibonding band. In Ref.liB it 
was argued that the latter is true, since the chain and the antibonding band anticross at the Fermi level. In 
this case reasonable numerical estimates lead to the conclusion that the total—tunneling current will have the 
desired phase shift, unless the chain band gap is too small. The authors of Ref.Eil made an opposite assumption, 
namely that the signs of the order parameter in the chain and the antibonding band are different. This makcsthc 
gap disappear near the band crossing and thus fits with various experiments indicating vanishing minimal gap. 
However, it is more difficult to explain the Josephson experiments in this model (one would have to assume that 
both in the chain and in the bonding bands gaps are much smaller than in the antibonding band). 

While the interband sign reversal model offers an alternative explanation of the Josephson experiments, it 
appears to be the only one able to explain the neutron scattering result in a plausible manner. Currently there is 
a consensus in the neutron scattering community that in fully doped YBa2Cu307 a sharp peak in the imaginary 
part of magnetic susceptibility x" appears below T c at the plane wave vector q=( — , ?), which shows perfect 
antiferromagnetic correlations between the planes. The most obvious effect of the onset of superconductivity 
on magnetic susceptibility is the appearance of a peak in x' at fiw = 2A, which however is allowed only if the 
correspondent wave vectors connects the parts of the Fermi surface with the opposite signs of the order parameter. 
Another effect is the appearance of a peak in x" a t = A + , where is the position of the van Hove 
singularity. In the RPA approximation the two peaks enhance each other, thus leading to a strong effect. However, 
the only way to ensure that the peak has no ferromagnetic (between the planesX-component is to assume that 
there is no intraband sign reversal at q=(2-, but only interband sign reversal.E3 

It is worth noting that intraband sign reversaLdoes not exclude any angular anisotropy of the order parameter. 
Moreover, the specific model considered in Ref.U appeared to have rather strong angular variation of the gap 
value, with the symmetry close to (x 2 — y 2 ) 2 . There are numerous experiments pointing to large variation of the 
absolute value of the gap in YBa2Cu307. As long as these experiments do not distinguish between d— and strongly 
anisotropic s— wave states we do not discuss them here. A final note concerns isotope effect, which is known to 
be small at the optimal doping (the highest T c ), but increases rapidly when superconductivity is suppressed by 
oxygen reduction, Pr doping, or other means. Such a behavior is very typical for superconductivity of mixed 
phonon and non-phonon origin. However, both components should work constructively, and not destructively. 
For instance, in a regular electron-phonon isotropic superconductor spin fluctuations, if any, are pair-breaking, 
and not pairing. One possibility to rationalize the isotope affect is to say that the electronic component is non- 
magnetic (like acoustic plasmons), or that the major process is two-magnons exchange. Another possibility is to 
say that the phonon- and spin fluctuation-induced interaction are separated in the phase space. For instance, if 
(as it is usually assumed) spin fluctuations are operative near q=(^, f ), but phonons only at small q's, both can 
be pairing in the contest of the d x 2_ y i superconductivity. Alternatively, if spin fluctuation are operative only in 
the interband channel, as it follows from the neutron scattering experiments, and phonons only in the intraband 
channel (in other words, only even [gerade] phonons interact considerably with the electrons), the isotope effect 
trends follow naturally. 
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